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We present a minisuperspace analysis of a class of Lorentzian 
wormholes that evolves quantum mechanically in a background Fried- 
man Robertson Walker spacetime. The quantum mechanical wave- 
function for these wormholes is obtained by solving the Wheeler- 
DeWitt equation for Einstein gravity on this minisuperspace. The 
time-dependent expectation value of the wormhole throat radius is 
calculated to lowest order in an adiabatic expansion of the Wheeler- 
DeWitt hamiltonian. For a radiation dominated expansion, the radius 
is shown to relax asymptotically to obtain a value of order the Planck 
length while for a deSitter background, the radius is stationary but al- 
ways larger than the Planck length. These two cases are of particular 
relevance when considering wormholes in the early universe. 
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1 Introduction 



Wormholes are handles in the spacetime topology linking widely separated 
regions of the universe, or "bridges" joining two or more different spacetimes. 
Current interest in Lorentzian or Minkowski-signature wormholes has been 
motivated initially by attempts to understand nontrivial spacetime topolo- 
gies, topology changing processes, their physical consequences, and their role, 
if any, in quantum gravity. A survey of recent work on wormhole physics 
includes topics addressing the fundamental properties of traversable macro- 
scopic wormholes |2|] , construction of explicit wormhole solutions and cor- 
responding classical || and quantum || stability analysis, wormholes as time 
machines and causality violation || , wormholes in higher-derivative gravity 
|| , the phenomenon of gravitational squeezing of the vacuum as a natu- 
ral mechanism for producing the states of weak energy condition violating 
matter required to prop open the throats of semiclassical microscopic worm- 
holes [|7| , wormholes as representing one possible mode of fluctuation of the 
spacetime foam || , and wormholes as gravitational lenses ||. 

Recently, a very simple model of classical dynamic Minkowski signature 
wormholes connecting Friedman- Robert son- Walker (FRW) spacetimes was 
constructed with the aim of initiating an investigation of the possible cos- 
mological effects of such objects in the early universe |10| . In that work, 
an alternative solution to the horizon problem was proposed based on the 
possibility that the early universe did not necessarily inflate, but was popu- 
lated with a network of microscopic Planck scale wormholes joining otherwise 
causally disjoint regions of space. The essence of the idea being that these 
wormholes can allow two-way transmission of signals between spatially sepa- 
rated regions of spacetime and thus permit such regions to come into thermal 
contact. This wormhole "phenomenology" is based on very elementary con- 
siderations. Indeed, to be of cosmo logical interest, these wormholes need only 
have stayed open long enough for the radiation, initially located in causally 
disconnected regions, to traverse the "handle" . Within the class of wormholes 



treated in |T0[, a crude estimate for the time required for thermalization can 
be had considering the following purely illustrative example. First, assume 
an initial number density n(tpi) of Planck-sized wormholes of average radius 
Lp at the Planck time and let n(tpj) = jnpi, where 7 = const, represents a 
fraction of the Planck number density: npi = (Lp) -3 . Then, at some later 
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time t, n(t) ~ R 3 , so the wormhole number density is given by 




n{t)=in Pl U£L . (1) 



Define the volume-filling factor (the volume of space filled by the interior of 
wormhole mouths) by 

number density x wormhole volume n(t) ( l(t)\ 3 
Planck density x Planck volume npi \Lp J 

The instantaneous wormhole throat radius l(t) = aR, where R is the scale 
factor of the background spacetime and a is the wormhole radius in the limit 
of a static background. For a radiation dominated cosmology and with a 
matter distribution on the throat subject to a particular equation of state 
(P is the surface pressure density and a the surface energy density on the 
throat) P = -P(cr), one can show that the throat decouples from the expansion 
with / ~ R 2 , so that the filling factor becomes 

«/*(*) =7 f^V- ( 3 ) 



R 



pi 



When v fm (T) rs 1 , there is one wormhole per unit volume by this equilibrium 
time T, and the average particle will have traversed at least one wormhole. 
So, for example, an initial wormhole number density of only 7 m 10 -14 allows 
thermalization of distant spatial regions by the GUT time scale Tqut ~ 
10 _34 sec. Of course, different values of 7 will lead to different values for the 
thermalization time scale T. 

While the phenomenological application put forward in |TIJ is novel, it is 
far from constituting a proof as such, and much more work is clearly required 
to settle the issue. As a case in point, note that a constant 7 implies the 
total number of wormholes is constant. This need not be the case in general 
and we expect 7 = 7(t). This quantity is related to the initial spectrum 
of wormholes and a first principles calculation of this factor would provide 
significant insight into the structure of spacetime at the Planck scale. Of 



equal importance is the fact that the wormholes treated in [IU| were strictly 
classical and if such objects do indeed exist at or near the Planck scale, it is 
much more likely they will be quantum mechanical in nature. Moreover, one 
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would like to know how the wormholes evolve in time and how the specifics 
of their temporal development is affected by the background spacetime. It is 
this latter issue which shall be taken up in the present paper. 

The quantum mechanics of Minkowski signature wormholes in FRW back- 
grounds will be addressed within the framework of a minisuperspace analysis 
applied to a specific class of wormholes resulting from surgically modified 
FRW cosmologies. The wavefunction of the wormhole will be obtained as a 
solution of the Wheeler-DeWitt equation. Because no one knows how to solve 
the Wheeler-DeWitt equation on the full extended superspace of three dimen- 
sional metrics (and matter field configurations), it proves extremely useful to 
examine this equation on minisuperspace configurations, which are restricted 
configurations described by only a finite number of degrees-of-freedom, with 
all the other modes of the full superspace frozen out. This is a standard 
approximation employed in quantum cosmology and is easily adapted for 
quantizing localized configurations in spacetime. Such analyses have been ap- 
plied to the study of quantum wormholes in Minkowski, Nordstrom-Reissner 
and Schwarzschild-deSitter backgrounds j|, jlTf . While the Wheeler-DeWitt 
equation is time-independent, the expanding background provides a natural 
time parameter which an external observer can employ in order to character- 
ize change in the wormhole's configuration. For FRW spacetimes, this time 
parameter is simply the scale factor (or any invertible function of the scale 
factor) . 

This paper is organized as follows. In Section II we describe in detail the 
minisuperspace model we adopt for Lorentzian wormholes. Canonical quan- 
tization of the wormhole is carried out in Section III and the Wheeler-DeWitt 
equation solved in perturbation theory for the wormhole wavefunction. The 
expectation value of the wormhole throat radius is calculated and is found to 
be partially dragged along with the expansion when the universe is radiation 
dominated, although relaxing to attain a value of the order the Planck length 
for late times. By contrast, a quantum wormhole in a deSitter background 
maintains a constant throat radius that is always greater than the Planck 
length. Both these statements hold when there is no matter residing on the 
wormhole throat. We briefly comment on the effects of including matter 
in Section IV. Discussion and caveats are presented in Section V. A Green 
function needed for the solution of a certain nonlocal differential equation 
is calculated in the Appendix. Units with G = c = 1 are used throughout, 
while h = Lp, where Lp is the Planck length. 
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2 Minisuperspace Model for Wormholes 



The wormholes we will be considering result from surgically modified Friedman- 
Robert son- Walker spacetimes. We adopt this technique to make the subse- 
quent analysis tractable, however, we assume the qualitative features of the 
wormholes are independent of the details of the construction. To construct 
them, take two copies M 1) M 2 of a FRW spacetime with identical scale 
factors R(t) and spatial curvature constant k: 



ds 2 = -dt 2 + R 2 (t) 



-^-^ +r 2 (d6 2 + sin 2 9 d<j) 2 ) 
1 — kt 1 



(4) 



and remove from each an identical four-dimensional region of the form f2 12 = 
{r 12 < a}. The resulting spacetime contains two disjoint boundaries 90i j2 = 
{^1,2 — a } which are timelike hypersurfaces. An orientation preserving iden- 
tification dQi = dVt 2 = dM yields two FRW cosmologies connected by a 
wormhole whose throat is located on their mutual boundary dQ. This model 
is an extreme version of a wormhole for which the spacetime curvature in the 
throat is much greater than that in the regions surrounding the mouths. Here 
a = a(r), is a function describing the dynamics of the wormhole's throat. It 
is important to point out that the proper motion of the wormhole throat is 
in general independent, that is, decoupled, from the background cosmic flow. 
Nevertheless, at any given instant, the physical radius (see the explicit form 
of the throat's three- metric below in (9)) of the wormhole is given by the 
product aR. The wormhole is spherically symmetric and the boundary layer 
is just the world volume swept out by its throat. This procedure also leads to 
a wormhole connecting a single FRW spacetime to itself if one identifies the 
two background geometries, i.e., we have a FRW spacetime with a handle. 
In this case, the two regions f2i j2 can be separated by an arbitrarily large 
distance in an open (k = 0,-1) universe. For the classical case, one insures 
that the modified spacetime is itself a solution of the gravitational field equa- 
tions by effecting a proper matching of the metric across the boundary layer 
[ |T3| ] . These classical wormholes must necessarily involve stress-energy distri- 
butions that violate the weak energy condition in the vicinity of the throat. 
This feature is, however, not an artifact of the construction and holds in 
general [I], H. At the quantum mechanical level, the state of the wormhole is 
described by appropriate solutions of the Wheeler-DeWitt equation, either 
with, or without the inclusion of matter. 



6 



As a consequence of the above construction, the Ricci tensor is everywhere 
FRW except at the wormhole throat: 




WW = n Ari + + 4^) e(-i7), (5) 



where R^^,R^^ are the Ricci tensors for the two spacetimes M 1 and M 2 , 
respectively and A/C* = (/Cj 2 ^ — /Cj 1 '' 4 ) is the jump, or discontinuity, in the 
extrinsic curvature of the boundary layer (i.e., the throat) in transiting from 
Mi to M.2- Proper normal distance as measured from the throat is conve- 
niently parametrized in terms of the Gaussian normal coordinate 77. For the 
case at hand, reflection symmetry implies A/C* = 2JCf^ 1 and spherical sym- 
metry implies that /C* = diag(/C£, K, e e , K, e ). The gravitational action for the 
wormhole is given by two contributions, to whit 



Swormhole Sthroat ~\~ Sbackgroundi (6) 

where 

Sthroat 



and 

^background 



— [ d 4 x^R^-—[ d A x^RV. (8) 

.67T J Mo 16lV J Mi 



167T JM 2 " 167T J Mi 

We shall first calculate the throat action. The 3-dimensional metric hij 
induced on the throat world-volume is 

hij(x) = diag (-1, a 2 (T)R 2 (t), a 2 (r)R 2 (t) sin 2 0) , (9) 

where r denotes the proper time on the wormhole throat and t refers to the 
background cosmic time. To evaluate (7), we need the components of the 
extrinsic curvature tensor 

ICij =%V (j )ef. } , (10) 

where the e(j) constitute a set of three linearly independent tangent vectors 
defined along the intrinsic coordinates £ l parametrizing the world-volume 
hypersurface and n M is the outward unit normal {n^n a = 1). The covariant 
derivative intrinsic to the hypersurface is taken along the jth coordinate 
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direction. The throat proper time and the two angles provide a convenient 
set of intrinsic coordinates: £ l = (r,0,(f>). As the location of the throat 
(i.e., its embedding with respect to the background FRW geometry) is = 
(t, a, 9, 0), the tangent vectors are given by = dX^/d^\ Since n^e^ = 
for i = 1,2,3, the unit normal to the throat is given by 



n 



aR 1 



— Vl - kg 2 + a 2 R 2 , 0, . (11) 



yi-/ta 2 ' i2 
A straightforward calculation of (10) yields 

aR' 1 

^4, = = —f- 
and 

di? 2di?' (na)d 2 R 
K, = —p H — H , (13) 

yl - Ka 2 + (aR) 2 Vl - ^a 2 (1 - /ta 2 ) yl - /ta 2 + (di?) 2 

where d = da/dr and it!' = dR/dt. Note the appearance of the two time 
parameters r and t in (12) and (13). We can always cast the equations in 
terms of one or the other time parameter by means of the jacobian 

_ dt _ /l - Ka 2 + (di?) 2 \ 1/2 
6r ~d^~ { T^a~ 2 ) • { ' 

Physically, this amounts to a calibration of the clocks attached to the throat 
in terms of the comoving clocks, if we choose to express r in terms of t, or 
vice-versa, if instead we choose to express t as a function of r. We shall 
have occasion to use (14) below. From (7), (9), (12) and (13), we have 
Tr{K) = A/C; + 2A/C^ and 

f , r , 9 9 , i ( aR \ 3aR'(a 2 R 2 ) / , 

Sthroat = J dr[2(aaR 2 + a RR) smlT 1 \ —=== \ - — L =( 15) 

- 2aRyJl - Ka 2 + (di?) 2 ], 
where we have made use of the identity 

dR d . / dR \ di? 

- /ta 2 + (di?) 2 dr Sm \^l-Ka 2 ) ^1 - no 2 + (di?) 2 

Kad 2 i? 



4±#) 



1 - /ta 2 ),/l - na 2 + (di?) 5 
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together with an integration by parts. Because R = R'(dt/dr), we are able 
to express R in terms of R,R',a and a. 

The contribution to the action coming from the background (external to 
the throat) spacetime is 

f dv V 2 2 /2 

S background = 3 dt , =(R R + RR +K.R), (17) 

where M = (M 2 — f2 2 ) \J(Mi — f^). Note that the integration over r is 
formally divergent for the cases k = —1,0, corresponding to open universes. 
To avoid possible technical complications associated with such divergences, 
we consider for the moment the case when the background FRW spacetime is 
spatially bounded, though the point we will be making below is independent 
of the value of k. In this case, the integral over r is convergent and gives 
(remembering that due to the spherical excluded region, a < r < 1) 



f 1 drr 2 ix a n ^ 1 . , ,. . , . 

L T7T=7i = 4 + 2^"+ 2 sln W s /(0) ' (18) 



»i drr 2 

so that 

\2\ 1/2 

... - 3 I n.A I 1 ~ " T ' ' 

background 



= sjdr f(a) °* fl 2 ( R2R " + RR ' 2 + R ^ ( 19 ) 



employing the time calibration jacobian in (14). With (15) and (19), we can 
now easily extract the wormhole lagrangian from 

Swormhole J dT Lwormholei.Q'i ^ j -^5 ^ )j (20) 

and calculate the momentum conjugate to a, II = dL/dd : 

\VT^J VT^a 1 - a 2 + (dR) 2 K ' 

(21) 

where 

F K (R, R', R") = R 2 R" + RR' 2 + kR. (22) 

As it stands, the transcendental relation (21) in general cannot be inverted 
in closed form in order to yield a in terms of II and a. The 'culprit' of 
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this obstruction is the jacobian factor (dt/dr), which depends explicitely 
on the throat velocity a for all values of k. However, provided there exist 
backgrounds for which F K = 0, then inversion of this relation in closed form 
is possible, and we can solve for the throat velocity in terms of the throat 
momentum (valid for all k): 

■ VI - na 2 . , / IT aR' \ 

a = smh H ; , (23) 

R \2aR 2 2^/1 -na 2 ) 1 ; 

where the canonical momentum IT is that of the wormhole's throat. For- 
tunately, there are several physically interesting cases of background FRW 
cosmologies for which F K vanishes identically. Indeed, they are (i) Minkowski 
space: k — 0, R — 1 ||, (ii) radiation dominated expansion: k = 0,R(t) = 

t 1 / 2 , and (hi) DeSitter space: k = 0,R(t) = eV™^', where A is a cosmo- 
logical constant. For this latter case, one of course adds a cosmological con- 
stant term to the action in (8). We shall consider wormhole dynamics with 
p K =o _ o £ or rema j nc l er f this paper. We should point out that cases 
(ii) and (iii) are, of course, of particular relevance when treating quantum 
wormholes in the early universe. Moreover, as the curvature k is negligible 
during the early stages of expansion, taking k = is a good approximation 
during this epoch. 

From the lagrangian in (20) the classical Wheeler-DeWitt hamiltonian 
for the wormhole is obtained by the standard Legendre transform : 

H(U, a; R, R') = Ud - L(a, a, R, R') (24) 



-2a 2 RR'Jl + (di?) 2 sinh _1 (di2) 



+ 2dR'(a 2 R 2 ) + 2aR^l + (dR) 2 

9 , / n aR'\ ( 
-2a RR cosh — — + — — — + 



\2aR 2 2 J \2aR 2 2 

n 

2ZR 2 



, nn , , / n aR' 
+ 2a 2 RR' sinh + 



„ , / IT aR'\ 
+ 2aRcos^— + — j, 

where in the final equality we used (23) to eliminate a in terms of II. As 
is to be expected, the classical dynamics of the wormhole follows by setting 
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H (IT, a; R, R') = 0, the vanishing of the hamiltonian being a consequence of 
the reparametrization invariance of the gravitational action. Note that this 
hamiltonian depends explicitely on external cosmic time through the scale 
factor R and its first derivative R' . Moreover, the canonical momentum of 
the throat (21) is shifted by an amount directly proportional to time rate 
of change of the scale factor. The one degree of freedom in this (classical) 
minisuperspace is the throat radius function a. 



3 Quantum Dynamics of the Wormhole 

As we are interested in the quantum mechanical features of Minkowski- 
signature wormholes, let us proceed directly to quantize canonically the class 
of wormholes constructed above. Invoking the prescription 

n— «£, (25) 

promotes the classical constraint H = to a differential equation 

H (^-ih-^,a;R,R^(a,t) = 0, (26) 

where ty(a,t) is the time-dependent wavefunction of the wormhole. As 
pointed out above, the throat momentum IT is shifted relative to the sinh -1 
factor by an amount proportional to the time rate of change of the back- 
ground spacetime, R'. This term represents the coupling between the worm- 
hole throat and the background. This shift results in a complex Wheeler- 
DeWitt hamiltonian operator 



H = H r + iHi (27) 



where 



H r = [-(2a 2 RR')x cosh(a;) + (2a 2 RR 1 ) sinh(x) + (2aR) cosh(a;)] cos($8) 

+ (2a 2 RR') sinh(:r) y sin(y), 

is the real part and 

Hi = [(2a 2 RR')x sinh(a;) (2a 2 RR') cosh(:r) - (2aR) smh(x)\ sin(yfl29) 

+ (2a 2 RR') cosh(x) y cos(y) 
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L 2 ft 

is the imaginary part. Here, x = aR'/2 and the operator y = ^rf^^r, 
with Lp = h 1 / 2 the Planck length, and we have adopted the factor ordering 
prescription such that y always stands to the right of x. In the absolute 
static limit (e.g., case (i)), R' = so x = 0, and we obtain (after setting 
R — 1, without loss of generality) 

h - = 2ac ° s (f 9 <3o) 

Hi = 0, (31) 

that is, the hamiltonian is purely real. Solutions of the Wheeler-DeWitt 
equation with the hamiltonian in (30,31) provide the starting point for a 
time-independent quantum analysis of wormholes in a Minkowski background 

I- 

Although we have not been able to find exact solutions of (26) for R' ^ 
0, we can obtain approximate solutions for an adiabatic expansion of the 
hamiltonian. To this end, we expand H r + iHi to order R', and put H = 
+ iR' HP where 

= 2ai?cos(y), (32) 
H\ l) = -(3a 2 i?) sin(y) + (2a 2 i?) y cos(y). (33) 

We solve for the wormhole wavefunction to O(R') in perturbation theory: 

^(a,t) = ¥ 0) (a,t) +R'¥ 1) (a,t) J (34) 

where is the exact zeroth-order solution which satisfies H^°^^(a, t) = 0. 
To this order, the Wheeler-DeWitt equation implies a relation between the 
zeroth and first order parts of the wavefunction: 

* W ( a> t ) = -i ¥V (a, t) . (35) 

We make use of the fact that there is an infinite set of exact zeroth-order 
solutions 

tt<?(M) = C m e- ( - m+1 ^ a2R2 / L p, (36) 

which is an obvious generalization of the wavefunction employed in [|J . 
The integer index m = 0, 1, 2, • • •, negative values being discarded as the 
wavefunction is not normalizable for m < 0. Since 

y^ ) (a,t) = -( m + 1 -)n¥°\a,t) : (37) 
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inserting (36) into (35) leads to 

cos(y)*«M) = R m+1 ^(o)( M ). (38) 
To solve this equation, we assume a product ansatz 

9£\a,t)=g(a,t)*<$(a,t), (39) 
which when substituted into (38) and employing the identity [T5j 

cos(y)g(a)f(a) = [cos(y)^(a)][cos(y)/(a)] - [sin(y)#(a)] [sin(y)/(a)], (40) 
leads to the equation 

3z 

sm(y) g(a,t) = -—a. (41) 

This nonlocal inhomogeneous differential equation can be solved using Green 
function techniques. The details of the calculation of the Green function for 
the operator sin(y) are relegated to the Appendix. The result is that 
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9W 



roc 

, )/T , J o G(u-u')Vtidu' + Y,A m e- m7ru , (42) 

where u = a 2 R 2 /L P 2 > and the Green function G(u — u') is calculated 
explicitely in (A10,A11). 

To proceed, we take the A m = (i.e., vanishing homogeneous solution) 
and note that the Green function (A10,11) has the form of a "kink" or topo- 
logical soliton with center at u > 0: 

{—\hk, for u' > u 
+wf , for u < u 

While (42) is exact, the integral cannot be reduced to known elementary 
functions. Nevertheless, a reasonable approximation can be had by taking 

G{u-v!) ~ J|(6(u-u') -e(u'-u)), (44) 

which does incorporate the essential kink structure and asymptotic behavior 
of G. Doing so, we obtain 

g (a, t) = -iV2^a 3 (R(t)/L P ) 2 , (45) 
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after discarding an infinite constant |T3[ . 

The wormhole wavef unction, to the order we are working, is therefore 
given by 

V mn (a, t) = ^l(a, t) (l - iV2^a 3 R'(R/L P ) 2 ) , (46) 

where 



-(™+2-Wf?) 2 ,-(™+D-(ff) 2 



(47) 



and is manifestly complex, the imaginary part being induced by the time 
dependence of the background spacetime. The leading order factor is sim- 
ply a linear combination of the exact zeroth-order solutions so chosen that 
^mn(O) = 0, ( m / n ) ^e condition required for hermiticity of [|J. 

With the FRW wormhole wavefunction in hand, we proceed to calculate 
the normalization constant C mn and the expectation value of the wormhole 
throat radius. From ^ 

/ da\y mn (a,t)\ 2 = 1, (48) 
J o 

we find after an elementary integration that 

C-2 = W * + ^ 2 ) (49) 

irr3 d/2 

+ ^Ff^ (( 2 ™ + i r V2 + (2" + 1)~ 7/2 - 2(n + m + l)- ?/2 ) . 



8vr 2 i? 3 

In a similar fashion, we find for the matrix element of the throat radius 

<ty mn \aR\^ mn > = < aR> mn = R<a> mn =R daa\ty mn (a,1$b) 

Jo 

T f 1 4- 1 21 



l2m+l 2n+l n.+m+lJ 
7T 1 1 1 2 



v / 2m+T 1 V2n+1 V"+"H-1 J 



7T Z /t 



where the numerical factor f( m ,n) is defined by 

1 I _^ 2 1 1 i_ 1 2 

(2m+l) 4 "r" (2n+l) 4 (n+m+l) 4 J it) (2m+l) 7 /2 ~T ( 2n+ i)T/2 ( n+m+1 )7/2 



J(m,n) J-^ i 



+ _J 2 4 1_ , 1 



(2m+l) 1 (2n+l) (n+m+1) VSm+l V2n+1 

(51) 
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Thus, to lowest order, the average value of the throat radius is of the order 
the Planck length with corrections coming from the explicit time depen- 
dence of the background spacetime. From (50), we can see that for static 
backgrounds, we recover the result that the wormhole is quantum mechani- 
cally stabilized against collapse, with an average radius = Lp 0. However, 
when the background depends on time, the wavefunction becomes manifestly 
time dependent and the wormhole throat radius need no longer be stationary. 
Indeed, already at lowest order, we find that the radius is modified from its 
static value. Thus, for radiation dominated expansion (RDE), (-^-) 2 = ^2, 
while for a deSitter cosmology, {^) 2 = |A|/3. As f( m ,n) > 0, we see from 
(50) that the expansion augments the magnitude of the matrix element. It 
is noteworthy that for the radiation dominated case, the perturbation to the 
throat radius expectation value vanishes as t — > 00 so that the wormhole 
throat relaxes from initially larger values to assume its static value. For the 
deSitter example, the first order perturbation is actually time independent; 
in summary: 

<aR> t -> < aR > statlc « L P , RDE, (52) 
< aR > t > < aR > s t a tic deSitter. (53) 

These calculations indicate that the expectation value for the wormhole ra- 
dius is larger than its static value when the background is expanding, indicat- 
ing that the throat gets partially dragged open with the expansion. This is 
perhaps not too surprising as the wormhole is gravitationally coupled to the 
background. However, the extent of the dragging cannot be known a priori, 
and we clearly see that it depends crucially on the form of the background 
scale factor. 



4 Comoving Wormholes 

When matter is added to the throat, even approximate solutions of the 
Wheeler-DeWitt equation are difficult to obtain. The variant of the WKB 
approximation scheme described in [15], designed to handle non-quadratic 
hamiltonians, fails to be of use for our hamiltonian because the classical 
energy relation 

H (n, a; R, R') = E (54) 
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with H as given in (24), cannot be inverted in closed form so as to ex- 
press II = U(E, a; R, R'), the latter relation needed for calculating the WKB 
wavefunction. Nevertheless, we can make an exact statement, at least at the 
classical level, concerning the state of motion and type of matter needed such 
that the wormhole throat is comoving, that is, coupled to the Hubble flow. 

The most general stress-energy tensor which gives rise to two identical 
FRW spaces connected by a 5-layer wormhole is 

Tftx) = SZ(x) S(rj) + e(rj) + Q(- V ), (55) 

where = is one of the standard perfect fluid source terms leading 
to (4) and 

S£ = ]im f drjT^x), (56) 

is the surface stress energy on the wormhole throat. In components 

£} = diag(-<7,P,P) (57) 

where a, P denote the surface energy and pressure densities. The throat 
matter action and lagrangian are 

Smatter = - J d 3 x \[—ho = -4tt J dr a 2 R 2 a = J drL matter . (58) 

The Wheeler-DeWitt hamiltonian for the combined system of wormhole plus 
matter H tota i = H throat + H matter , where 

H ma tter = -^matter = 4:71a 2 R 2 a. (59) 

In order that the wormhole be comoving with the background expansion, 
a — 0, since the radius l(t) = aR(t) is then simply proportional to the scale 
factor. This is also reflected in the fact that dt/dr = 1 if and only if a is 
constant, i.e., throat proper time is identical to cosmic time when there is no 
relative motion between the throat and background. For a = the classical 
hamiltonian constraint, obtained after summing (24) and (59), reduces to 

2aR + 4ira 2 R 2 (T = 0, (60) 

that is, the surface energy density on the throat must be 

a = --^—<0 (61) 
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which is negative definite and scales as R~ l in time. According to this, the 
wormhole mouth will expand just as rapidly as the Hubble flow and for all 
scale factors. However, in order for this motion to be realized, the matter on 
the throat must be "adjusted" to have the surface energy density specified 
in (61) as well as obey the equation of state P = — | flO| . We suspect this 
will also hold at the quantum level, that is, a quantum comoving wormhole, 
one for which < aR >= constant x R, will require nontrivial stress-energy 
located at the throat. We can turn this argument around and conclude 
that in general, the wormhole will not be comoving, but will exhibit motion 
relative to the background expansion. 

This result, that a classical comoving wormhole entails weak energy con- 
dition (WEC) violating matter should be compared to the results of Ro- 
man |17|, who analyzed a class of Morris-Thorne type metrics represent- 
ing Lorentzian wormholes embedded in a deSitter inflationary background. 
Those classical wormholes were shown to inflate in step with the background 
and require a non trivial stress energy tensor violating the WEC. 



5 Discussion 

In this paper we have investigated a simple model of Minkowski-signature 
wormholes attached to FRW spacetimes. The quantum gravitational dy- 
namics of these wormholes is reduced to the quantum mechanics of a single 
variable, namely, the throat radius, by solving the Wheeler-DeWitt equation 
over minisuperspace. When quantizing wormholes in static backgrounds via 



the Wheeler-DeWitt prescription 0, [II], one comes up against the problem 



of time. That is, the resulting wavefunction is time-independent, despite the 
fact that the wormhole may well evolve classically. This makes it awkward 
to discuss the quantum stability since the calculation of transitions between 
allowed quantum states of the wormhole lies beyond the scope of the formal- 
ism. However, for time-dependent backgrounds one can treat the evolution 
of quantum wormholes in an obvious and natural way. For the case of FRW 
backgrounds, the temporal parameter is provided by the scale factor. More- 
over, FRW cosmologies are the relevant spacetimes in which to analyze early 
universe wormholes. 

We have calculated the average wormhole radius for both radiation dom- 
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mated and deSitter inflationary scale factors. It is noteworthy that even in 
the absence of matter at the throat, the wormholes are stabilized quantum 
mechanically against both collapse and blowing-up. For radiation domina- 
tion, the wormhole radius decreases monotonically from intially large values 
to attain a minimum of order the Planck length. For the inflationary case, 
the lowest-order calculation reveals the radius to be held fixed at a value 
always larger than the Planck length, by an amount proportional to the ab- 
solute value of the cosmological constant (vacuum energy density). These 
two cases indicate that the wormhole throat gets dragged open by the back- 
ground expansion, although there is relative motion, or "slippage" , between 
the wormhole and the ambient Hubble flow. This result is encouraging as it 
weighs in favor of the ideas set forth in [|TtJ . 



A limitation of our calculation is the fact that the wormholes treated 
here have zero-length handles. Of course, we by no means claim that all 
(Planck-scale) wormholes would be of this form, but we believe that the 
qualitative aspects of their evolution should not be overly sensitive to this 
trait. Nevertheless, it would be interesting to check this by quantizing finite 
handle wormholes as described, for example, by Morris-Thorne type metrics. 

Perhaps a more severe limitation of our calculation is due to the fact 
that we have not employed the full superspace description for the wormhole 
dynamics. The truncation from an infinite number of degrees of freedom 
down to a single degree of freedom may appear drastic, but it is the only 
viable option if we wish to be able to make a definite calculation. Working 
in the full superspace is obviously preferable, but is currently an intractable 
problem. 

Finally, when formulating a time-dependent quantum mechanics of worm- 
holes via the Wheeler-DeWitt approach, as we have done here, we find the 
hamiltonian operator picks up many new terms not present in the static 
background limit, e.g., compare (28,29) to (30,31). These additional terms 
are functions of the expansion velocity R' and complicate the task of solving 
for the wavefunction. Nevertheless, we have been able to calculate the wave- 
function perturbatively employing R' as a small parameter. Naturally, one 
would like to get at information lying beyond the restraints of perturbation 
theory. To this end, we would like to remark that the wave equation (26) may 
lend itself to approximate, but non-perturbative, solutions based on methods 



borrowed from the calculus of finite differences 181. This follows from the 



observation that the nonlocal operators appearing there can be expressed in 
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terms of finite translations acting on minisuperspace. These points will be 
explored elsewhere. 
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A Appendix 

We wish to obtain an explicit representation for G, defined as the inverse, 
over C°° function space, of the nonlocal differential operator sin(J^), that is, 

nn(J^G(u,u')=8(u-u'), (Al) 

where 

. ( d\ _ d Id 3 Id 5 

Sm [du) ~^~3W + 5!^~"' [A2) 

As the operator (A2) is linear, it proves useful to introduce the Fourier trans- 
form G of G via 

G{u, u) = r * G( p ) e HP(«-«') > {M) 

J-oo V 

whereby using the identity 

sin (j^j e- ipu = sin{-ip)e- ipu , {AA) 

we have that the transform of the Green function is 

G{ P ) = {A5) 
smh(p) 

The problem of computing G thus reduces to carrying out the integration in 
(A3) with G as given in (A5). There are two cases to consider: 
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A.l case (1): A = (u - u') > 

We go to the complex p-plane and consider the following contour built up 
from closed intervals and semicircular arcs centered about the origin: T^s — 
[-R, —5] + 75 + [5, R] + jr, where 75 is a semicircle of radius n > 5 > 0, 
and 7/j a semicircle of radius R with (M + 1)tt > R > Mir (M » 1), both 
located in the lower half-plane, as indicated in Fig. A. The inequalites are 
imposed so that the contour avoids hitting any of the poles. Since sinh(x + 
iy) = sinh(x) cos(y) + i cosh(x) sin(y), the integrand in (A3) has simple poles 
located on the imaginary axis: z = ±imir for m = 0, 1, 2, • • •. Thus, by the 



residue calculus [|16|1 , we can calculate the integral taken around the complete 



contour in the clockwise sense 



r R>s V2¥ smh(z) ^ 



The Fourier integral representation for G (A3) results after taking the two 
limits 5 — > and R — > 00. Denoting by Ir the contribution to (A6) coming 
from the large semicircular path 7^, it is straightforward to show that 

r w d9 g-A/?sin6» 

0<\I R \<R - ? =- -j- . — >0, (A7) 

Jo V27T I cosh (RcosO) — cos 2 (Rsm9)\ 1 ' 2 

as R — > 00 for every 9 G [0,7r]. On the other hand, the contribution Is to 
(A6) coming from the small semicircle 75 centered about the origin gives 



Putting together (A3) and (A6-8), we have that 

v 7 ^ E (-e^ A ) m = / 4=^TT + i im n Ja + i im Ir > ^ 
^ 1 J -00 \/2tv sinn p fl^oo 



m =i --^ v 27rsinh(p) 

so that 

where in the limit M — > 00 , i? — ► 00 and 5-^0 the contour maps to 
the real line joined to an infinite semicircle in the lower half-plane. 
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A.2 case (2): A = (u - u') < 



This time, we close the contour in the upper half-plane and in the counter- 
clockwise sense, i.e., T R ^ = [-R, —8} + 75 + [5, R] + 7^, where now 7,5 and 7^ 
are semicircular arcs in the upper half-plane. The steps and limits involved 
in computing G are entirely analogous to those of the previous case, with the 
final result that 



The Green function (A10,11) is bounded for |A| — > 00 and G(0) = 0. 
From (Al), the general solution of sm(jjj g(u) = f(u) is therefore 





gill) = ghomogeneousiu) + / d,U G{u - U ) /(it'), 



(A12) 



where the homogeneous part of the solution is given by 




0413) 
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Figure Caption Fig. A. The contour for case 
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